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Abstract 

We evaluate the multiplicity of hadrons in the e'''e~-annihilation at a given 
thrust T in the modified leading-log approximation, including 0(y^a7) corrections. 
The calculation is done at a large value of r = 1 — T by the use of the factorisation 
which takes place in the one-particle- inclusive cross section at a given r. At a small 
r, a different type of factorisation takes place, which also enable us to evaluate the 
multiplicity. Two approaches are compared numerically. Measuring this quantity 
near r = 1/3, we can determine the multiplicity ratio between a gluon-jet and a 
quark-jet. 



1 Introduction 



The Modified Leading- Log Approximation (MLLA) |T|[0 with the Local Parton- 
Hadron Duahty (LPHD)[^ has been successful in describing the majority of the e~^e~- 
annihilation events (see, for example, The shower Monte Carlo programmes 

based on it [0 are widely used for the simulation of experiments. Nevertheless, our ex- 
perimental knowledge of the multiple hadroproduction in gluon jets is still poor, or even 
apparently in contradiction with the theoretical predictions. The contradiction has been 
felt particularly in the multiplicity ratio between a gluon-jet and a quark-jet. 

In earlier time, some argued that the multiplicity ratio would be identical to the ratio 
between the colour charges of a gluon and a quark 0. Nowadays, this expectation is 
justified in the following way: According to LPHD, the average number (multiplicity) 
of hadrons produced in a hard process is proportional to the number of gluons emitted 
perturbatively from the hard partons. Because of the infrared singularity of the emission 
amplitude, the majority of the gluons are soft (much less energetic than the hard parton). 
At the leading order in MLLA, the strength of the soft-gluon radiation is proportional 
to the colour charge of the hard parton {Ca = 3 for a gluon, Cp = 4/3 for a quark) 
with which the jet is associated. Thus the multiplicity ratio between a gluon-jet and a 
quark-jet is expected to be close to Ca/Cf = 9/4 at high energies. 

The ratio has been measured using three-jet events of the e"'"e~-annihilation, iden- 
tifying one of the jets as a gluon-jet. The ratio was found much lower than the above 
expectation Ca/Cf- For example, OPAL collaboration gave 



(including neutral particles) and 

('^CH)g-jot 



1.267 ± 0.043(stat.) ± 0.055(syst.) (1) 



1.326 ± 0.054(stat.) ± 0.073(syst.) (2) 



(^CH)q-jot 

(charged particles only). (For earlier experiments, see [p^ .) 

The next-to- leading order correction in MLLA {0{y/al) = 0(1/ ^\n{W'^ / Aq^^j^)), 
W: the total energy) reduces the leading order prediction r = 9/4 by about 10 percent 



TT| . The correction is far too small to explain the discrepancy from the experimental 
observation. (The next-to-next order correction does not change the situation [0].) 

It should be noted, however, that the multiplicity in a gluon-jet {n)g_j^^, which was 
calculated in MLLA, is not necessarily identical to the quantity measured in the experi- 
ments. In the theoretical calculation, it is defined as a half of the total multiplicity from 
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the two hard gluons created by a gauge invariant gluon source (for example, F^^, where 
F^y is the field strength of gluon). 

In the experiments, on the other hand, one selects three-jet events (defined one way 
or another) and compares the number of hadrons in each angular region. At current 
energies, however, a hard parton does not necessarily produce a well-coUimated jet, and 
it may cause biases, depending on how the measurement is made: If one does not impose 
the condition that each jet be isolated (i.e. well-collimated) , the contribution from each 
hard parton may mix up one anotheiQ. Thus the difference in the multiplicity due to 
the difference in the colour charge of the respective hard parton may be reduced in the 
comparison of the multiplicity in the respective angular region. If one selects the events 
with well-collimated jets, on the other hand, one is essentially comparing the multiplicity 
in the direction of the respective hard parton (in a narrow angular cone around the jet 
axis). Now the stronger radiation of gluons from a hard gluon, which is responsible for 
a larger multiplicity in a gluon jet, also causes a wider angular distribution of hadrons, 
owing to the larger recoil from the multiple gluon emission from the hard gluon @] (see 
the angular distribution of hadrons measured in 0). The multiplicity counted in the 
direction of the hard gluon, therefore, is not as large as Ca/Cf times the multiplicity 
counted in the direction of the hard quark or antiquark (see also 0). 

Though the magnitude of the biases is in general hard to estimate theoretically, we 
may examine it using the MLLA-based shower Monte Carlo simulation of the exper- 
iments. The simulations done by the experimentalists (in p and |jlO|) in fact gave 
reduced multiplicity ratio. It suggests, at least, that the observed smaller ratio does not 
necessarily contradict QCD itself. 

In any case, a bias free measurement of the multiplicity ratio is not easy, because a 
gauge-invariant two-gluon source is hard to prepare^. 

In this article, we analyse the multiplicity at a given value of thrust, where thrust T 
is defined, in the cm. frame, by 

T = max|^^^| {ff = l) (3) 

(the direction of the three- vector n is chosen to maximise the rhs). T = 1 would imply 

that all the particles are parallel (or antiparallel) to n. 

^ Strictly speaking, we cannot tell, in a gauge invariant way, which particle is emitted from which 
hard parton, except in the high-energy limit where only coUinear particles are produced. 
^The process 77 gg may be used for this purpose 
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If T is far from one, it implies that a hard gluon is emitted. The emitted hard gluon 
causes an increase of the multiphcity. At high energies, where most of the events with 
large r = 1 — T(< 1/3) are well-collimated 3-jet events, the multiplicity increase can be 
identified with the contribution from the multiplicity in a gluon-jet. In this way, we can 
avoid the experimental bias caused by the angular restriction, and we can compare the 
theoretical calculations with the experimental data without ambiguities. 

In Sect. 2, we evaluate the multiplicity at a large r at the next-to-leading order of 
MLLA. At this order, a hard gluon adds the multiplicity independently of the quark- and 
the antiquark-jets. The multiplicity increase is thus close to Ca/{^Cf) times the multi- 
plicity from the gg-pair, which is identical to the multiplicity (n)e+e- (without specifying 
the thrust value) at the same cm. energy as the gg-pair. 

The approximation, however, is not accurate if the energy is not high enough or r 
is not large enough. In either case, the gluon jet is not well-isolated, and higher order 
corrections cause substantial modifications. 

In Sect. 3, we evaluate the multiplicity at a small r, where no hard gluons are emitted. 
At the leading order, the multiple soft-gluon emission from the gg-pair determines the 
value of 1) 0||]15[. Because the soft gluons are emitted independently from the 
hard gg-pair, the one-particle-inclusive amplitude, which gives the multiplicity, factorises 
from dcr/dr, and the resummation to all orders can easily be done. 

In this approximation, we do not separately evaluate the multiplicity from a gluon-jet 
and a quark-jet. When we apply the approximation to a larger value of r, however, we 
expect two distinctive contributions, one identified as the multiplicity from the gluon-jet 
and the other from the quark- and antiquark-jets. The approximation is not justified 
near r = 1/3, because it assumes that the gluons are emitted independently. The error, 
however, would not be large in the intermediate r region. We therefore expect that 
the prediction in this region is close to the next-to-leading order prediction obtained in 
Sect. 2, if the energy is high enough. 

We evaluate the multiplicity in the small r region making use of the detailed ex- 



perimental data of the thrust distribution at 58GeV [20|. The multiplicity evaluated 



in this approximation in the intermediate r region (r ~ 0.1-0.2) turns out to be lower 
than the prediction in Sect. 2. This is because the multiplicity from the hard gluon is 
not well-separated from the multiplicity from the gg-pair. It suggests that we may have 
to go to higher energies in order to be able to determine the multiplicity ratio from the 
multiplicity measured at r in the intermediate region. 
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2 Multiplicity at large r (Next-to-Leading Order) 



The one-particle-inclusive (IPI) cross section counts the number of particles at a given 
momentum. When we integrate it over the momentum, and divide it by the cross section 
(which counts the number of the events), we obtain the number of particles per event. 
The multiphcity at given t — 1 — T is obtained by dividing the integrated IPI cross 
section at given r by the cross section dcr/dr. 

When T is large (but < 1/3), a hard gluon is emitted. At high energies, the multiphcity 
is given as the sum of the two distinctive contributions: 

(a) The registered particle originates from the hard g^-pair; or 

(b) it originates from the hard gluon. 

In either case, the registered particle is mostly produced near the direction of the respec- 
tive hard parton. 

At the leading order, the multiplicity at a large r is simply 1 + Ca/{'2Cp) times the 
total multiplicity (n)e+e- at the same energy. The additional multiplicity proportional 
to CaI{'2'Cf) is due to (b) (the multiphcity in a gluon jet). 

We first evaluate the next-to-lcading order corrections {0{'^/al)) to (a). Let Pi (P2) be 
the momentum of the hard quark (antiquark), P3 the momentum of the hard gluon, and 
k the momentum of the gluon to which the registered hadron belongs (Fig. la). Suppose 
that the hard gluon is emitted on the quark side (in the cm. frame). 

At the next-to-leading order, only the soft and/or collincar emission contributes to 
the multiplicity. We can therefore assume 2Pi ■ P3 ^ m^, where = 2Pi • k if the 
A;-gluon is emitted on the quark side, and = 2P2 ■ /c if it is emitted on the antiquark 
side. The value of thrust is determined by 2Pi ■ P3: 

2Pi-P3>m2 . 

Owing to the strong inequality, the cross section for the emission of the two gluons 
factorises into the cross section for the hard-gluon emission (dcr/dr) times soft- and/or 
coUinear-gluon emission amplitude. Therefore, when we divide the IPI cross section 
(illustrated in Fig. la) by dcr/dr, we obtain 

.rW'VSdm^ /•! 2Cf f I , Xi\ . 2x,./ 2x 

< >gg = / — - / dxi 1 + OLs[xxm )Mg{xim ) , (5) 

Jo m? J'm?/W'^ TT \Xx 2 / 

= 2Pi • P2 , 
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where Mg{Q'^) is the multiphcity from a gluon emitted at the transverse momentum Q, 
and the integration kernel represents the probabihty of the emission of a gluon at the 
longitudinal momentum fraction xi (cf. (10) below). The invariant mass of the gg-pair 
W depends on the momentum of the hard gluon P3. Later, we shall integrate the rhs of 
(5) over P3 with the weight of the hard-gluon emission probability. 

in the 



At the next-to-leading order, MLLA gives the multiplicity function Mg{Q 



form §|T6| 



M,(Q2) = C (in(QVA'))'' exp \2^J\n{Q^/KlaD 



(6) 



where 



7o 



7i 



ho 




Ca 



127r 



{Nf. the number of the active quark flavours). The normalisation constant C in (6) is not 
determined in the perturbation theory, and has to be fixed by comparing the predictions 
with the experimental data. The effective QCD coupling asiQ"^) becomes large as Q 
approaches the QCD mass scale Aqcd, and the interaction becomes nonperturbative. 
The emission of the strongly interacting gluons at low Q is responsible for the bound 
state formation, and does not contribute to the multiplicity. Mgi^Q"^), which represents 
the number of perturbatively emitted gluons, therefore, should be set zero for < Ql 
in the integration on the rhs of (5), where the infrared cutoff Qo is of the order of 
Aqcd- Qo is an unknown parameter from the point of view of the perturbation theory. 
The dependence on Qo, however, is absorbed in the normalisation constant C of the 
multiplicity function Mg-. i.e. the ambiguity disappears in the physically observable 
multiplicities owing to the renormalisation of Mg. 

The upper bound of on the rhs of (5) (m^ < tW'^/2) is an ambiguous quantity in 
our approximation. We have assumed at the accuracy of the next-to-leading order that 
tW"^ ~ 2Pi ■ P3 ^ m^. It implies that the upper bound of can be anywhere below 
tW"^ (if not much smaller than it), and that the dependence on the upper bound is indeed 
at the next-to-next order. This fact can easily be checked by the explicit evaluation of 
the rhs of (5). 
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The integration on the rhs of (5) gives at the next-to-leading order 




?,^\uWyKlcD ^4 Ca \ Ca 



7o ^/ /^C'f 




Let us next evaluate the other contribution (b) , in which the registered hadron origi- 
nates from the hard gluon. In this case, the multiplicity function Mg is attached to the 
hard-gluon line (momentum P3) (Fig. lb). Other gluons may be emitted from the ^fg-pair, 
but they have to be soft and/or collinear gluons at the next-to-leading order in MLLA. 
For the unregistered soft and/or collinear emission, we should note the following: 
• The soft and/or collinear emission does not change the value of thrust. 
■ The soft and/or collinear emission amplitude factorises. 

It implies that the singularity of the real emission amplitude is cancelled by the singularity 
of the virtual correction. Namely, the same cancellation that occurred in the case of the 
IPI cross section without specifying the thrust value takes place also in this case. We can 
therefore neglect the emission of the unregistered gluons at the next-to-leading order. 

Suppose that the event consists in only three (on-shcll) hard particles: a quark at a 
momentum Pi, an antiquark at P2, and a gluon at P3. The 0{as) cross section for the 
single gluon emission is given by 

^ I 2Pi-P3 2P2-P3 2(2Pi-P2) 1 d^Pg 

da - aoSnCa. \i2P,.P,)W^ + {2P, ■ Ps)W^ + (2^ • P3) (2P2 • P3) / (2^ ' 
= (Pi + P2 + P3)' , 

where ctq is the lowest order cross section without any gluons emitted (7* — > qq). 

We assume for the moment that the antiquark is more energetic than the quark (in 
the cm. frame) (P2 > P{*). We take the direction of the antiquark as the negative ^-axis. 

The direction of the thrust axis (n in (3)) is identical to the direction of the most 
energetic parton. We therefore have to distinguish the two cases. 

(i) If the antiquark is more energetic than the gluon (P2 > P3), the 2;-direction is 
identical to the thrust axis. We then have 

2Pi • P3 = rW^ , 
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r2 

1 



2P2 ■ P3 = x{l - t)W' 
2Pi ■ P2 ^ {1 - x){l - t)W^ . (9) 

where x = P^/{Pi + Ps)^ {P^ = (i^" + P-)/V2). 

The differential cross section for the process 7* — > qqg is in this case 

da Cpa, \ 2Pi ■ P3 , 2P2 ■ P3 , 2(2Pi ■ Ps)!^^ | 



drdx 27r i2P2-P3 2Pi • P3 (2Pi ■ P3)(2P2 • P3) 

C^a, f T x(l -r) 2(1 -x)] 
= ^0-^^ — 1 r ^ 



27r [ — r) r tx \ 

^ao^Mr^x). (10) 

The condition Pa^ > Pg" imphes 2Pi • P3 < 2Pi • P2: namely, 

.<1^, (11) 



while P2" > Pi" imphes 2Pi • P3 < 2P2 • P3, or 

T 



1 -r 



< X . (12) 



(ii) If P2 < P3 , the direction of the hard gluon is identical to the thrust axis. We 
then have 

2Pi • P3 = ^~/~^ W' 
1 — X 

2P2 • P3 = -^tW^ , 

1 — X 

2Pi . P2 = rW^ , (13) 
da do Cpas r 2Pi ■ P3 2P2 ■ P3 2(2Pi ■ P2)W^ \ 



and 



drdx l-x 27r i2P2-P3 2Pi ■ P3 {2P^ ■ P^){2P2 ■ P^ 
(To Cpas { XT ^ l-x-T _^ 2 xf 



1 — X 27: [1 — X — T XT x{l — X — t) j 

^ao^Mr,x) (14) 

The condition P^ < P^ imphes 2Pi • P3 > 2Pi • P2: namely, 

X < ^ , (15) 

i — T 
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while P2 > Pi now implies 

\^<x (16) 
1 + r 

Finally, the case P° > gives the same cross section (now the direction of the quark 
is taken as the negative z-direction), and thus simply doubles the contribution. 

When a gluon is emitted at a large angle, the angular ordering of the succeeding soft- 
gluon emission is not exact at the next-to-leading order. Accordingly, the multiplicity 
from the Pa-gluon is not simply identical to Mg{Q'^), where is the transverse compo- 
nent square of P3. Instead, it is given by Mg{2Pi ■ P3) + Mg{2P2 ■ P3) - Mg{2Pi ■ Ps)^ 
The multiplicity from the hard gluon {n)g at a given thrust value (0 ^ r < 1/3) is then 
obtained by averaging it with the cross sections (10) and (14) over the range of x given 
by (11), (12) and (15), (16) respectively: 



<n>g 



J dx {Mg{2P, ■ P3) + M,(2P2 • P3) - M,(2Pi ■ P2)} ^^^^^^^^^^ 

1 1-2t 

j^' dx {Mg{TW') + Mg (x{l - r)W^) - Mg ((1 -x){l- r)W^) } 
x/i(r,a;) 



x/2(t,i) , (17) 



where 



A(r) = Pi(r) + P2(r) , 

1-2t 

Pi(r) = dx/i(r,x) 

1 — T 

_ + 8r - 3 - 2r + 2 , 1 - 2r 
~ 2r(l - r) r(l - r) 

F2{t)= f~dxMT,x) 



In- 





1-2t 




1-T 




I 

— T 






1 




1 


— r 



(l-r)2 1 (l-r)2 1 , 1 + T 
In ^- ^ - - In ^ — ^ ^ + 2 In ■ 



2(1 -r) 



/i and /2 are defined in (10) and (14). 

In order for this expression to be correct, P3 has to be hard 
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Let us now go back to {n)qq. Averaging the rhs of (7) over the hard-gluon momentum, 
we obtain 

Ca \ 3^1nH^VA|cD ^ Ca\ Ca JJjJ da/dr 



Ca \ S^lnWyAlcD ^ Ca\ Ca J J ] A{t) 
X I dxMg ((1 -x){l- t)W'') h{r, x) + dxMg{TW'')h{T, x) 



(19) 



The multiphcity at a given r is the sum of the two contributions: 

< n >T=< n >qq + < n >g (20) 
The thrust-dependent multiphcity ratio R{t) is defined by 

where 

1 . ^ 

< n >q= -<n>qq . (22) 

The expressions (17) and (19) can be expanded around r — 1/3. At 0(^/a7), we 
obtain 

< - >. = 1 1 + . "° ^1 M, (^] , (23) 



^ 1 1 + ^0 f 1 _ ^ 

Ca \ S^lnWyAlcD ^ Ca\ Ca 



r2' 



xU. mI^), (24) 



5i(t) = j]f dxh{T,x)\n^ 



l-T 

1-2t 

l-T 



3tx 



+ /;-' dx/,(r,x)ln^^^-^, (25) 



l-T 
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B- 



1-2t 

,(t)= /;-^ dx/i(T,x)ln(3(l-T)(l-x)) 

1-2t 

+ /^"^ da;/2(T,a;)ln(3T) 



The ratio R{t) is thus 



X J 1 + 



70 (Sf 



We can easily perform the integral on the rhs of (25) and (26) to obtain 
Si(r) = ln(3T)Fi(T) - F^ir) + ^F^{t) + InS F2(t) 

B2{t) = In (3(1 - r)) F,{t) - ^F,{t) + ln(3r)F2(r) , 



where 



F3{r)= ll'^ dx/i(T,x)ln- 

1— r 



^ {i(l - 2r)(3 + 2r) in 1^ + I(x - 4) in ^ 
+i(3r - 1)(7 + X) - \(r- - 2r + 2) m il^ln i^' 



1 — -i 

^4(t)= /_7 dx/2(T,x)ln- 

1 + T 



1 -r [ V(l Vl + ryj I'^Vl-ry Vl + r 
+ ^+^'lnn .llnliim 1 + l-2r (l-2r)(l + r) 

^"(1-^)^" (137)^- 27(1^ ^"TT7^" (i-r)^ 

/■t^ 1 
^^^'^^ " /^^ dx/i(T, x) In 



l-T 



2 + 3t, t 3-6r ,l-2r (l-3r)(5 + 3r) 
In- rln ^ ' 



l-T l-T T(l-r) l-T 2t(1-t) 
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Fair) 



l-2r 

1-T 



1-r 
1 + r 



dx/2(r, x) In 



1 



1 



r 



Sp( 



r 



-2 In- -In 

1 — T-^ 



— Sp 

1 + r 
2(1 - 7 



(l-x)2 

1 - T 



1, 2 1 
+ -In^ 

2 1 + r 



+ In r In 



r(l + r) 
(l-r)2 



2(1 



r 



r(l 



Sp 



1 - 2r 



Sp 



1-2t 



i^7(r)= /J;^ dx/2(r,x)ln. 



i+. 1 - a; - r 

-2 ^3 



'■'rT7"'Tr37)J + MSp(— )-Sp(i-r)} 



2(1 -r) 1 + r (1-r 

+ In ^ I In ^ + 'l!,,"^ ^ . I Sp ( — ^ 1 - Sp 



2r2 2(1 - r) r(l - r) [ V (1 " ^) V V 1 + r 

l + r2 (i_2r)(l + r 



Sp(a;) is the Spence function defined by 



^ dr 1 

Sp(x) = / — In^— . (31) 



/ox 1 — X 

In the hmit r 1/3, we find that 5i(r)/A(r) ^0, (z = 1, 2), and 

i?(l/3) = r , (32) 

where r is the conventional multiphcity ratio between a quark-jet and a gluon-jet calcu- 
lated at the next-to- leading order [p!T[| : 



We show in Fig. 2 the prediction for {n)r divided by (n)e+e- at Vr=58GeV. {n)e+e~ is 
given by 

/•! ^ 2Cp 
< n >e+e- = / — - / _ _ dxi 



1 


- 1 + 


Xi\ 




— a 






2 ) 




AT; 






Ca 
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In Table 1, the prediction at 58GeV is compared numerically with that at 91GeV. We 
find that the differences are negligibly small. 

In Fig. 2, we also show -R(t) at W=58GeV. It reduces appreciably as r decreases, while 
change in {n)r is more moderate. As the emitted gluon becomes harder (at larger r), it 
carries away a portion of energy from the gg-pair, and the increase in {n)g is partially 
compensated by the decrease in {n)gq in their sum {n)^ (recoil effect). 

The multiplicity ratio R{t) is in general not a directly measurable quantity. What 
we measure without ambiguity is not {n)g or {n)g itself, but their sum {n)^-. However, 
-R(l/3) = r can be determined experimentally in the following way: First, we measure 
{n)^- at r = 1/3. (In practice, because of the lack of statistics near r = 1/3, we measure 
{n)r for r < 1/3, and extrapolate it to r = 1/3, assuming a smooth dependence on r.) 
Now we note that {n)qq at r = 1/3 is identical to {n)e+e- at the cm. energy W/V3, 
which is the cm. energy of the qq pair (compare (7) with (34)). We can therefore use the 
experimentally measured {n)f.+e-- at the cm. energy W/V3 for {n)gq at r = 1/3. Then 
{n)g is defined as (?t,)t-=i/3 minus {n)gg, and R{l/3) is obtained by (21). 

It should be noted that the above result (23) and (24) cannot be applied to r > 1/3, 
because the 0{as) cross section vanishes for r > 1/3. 



3 Small T region 

So far, we have analysed the large r region, assuming that the event consists in 
three well-separated jets (i.e. the invariant mass of each jet is much smaller than W). 



When r is small, multiple soft-gluon emission determines the value of r and it is 
no more possible to isolate a single gluon jet. In this case, however, the soft gluons can 
be regarded as emitted independently from the initial hard gg-pair, and the soft-gluon 
emission amplitude factorises (as is the case with soft-photon emission in QED). Making 
use of the factorisation, we can evaluate the multiplicity at a given r without dividing it 
into the contributions from isolated jets. 

Let us first remind ourselves that if all the emitted gluons are soft r can be expressed 



as the sum of the contribution of each gluon |ll4[| : 



r ?a ^ , = min{2Pi ■ ki, 2P2 ■ h} , (35) 

i 

where ki is the momentum of each soft gluon. The emission amplitude of n soft gluons is 
factorised into the n identical amplitudes of one-gluon emission (Eikonal Approximation), 
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and we obtain the thrust distribution F{t) in the form of a series 
F(r) = 

a dr 

= A{w', Ql) [rirW) + ^ £ dnr{T,w')r ((i - t,)w') + ■ ■ 

(r^O) , (36) 
where T^rW'^) is the factoriscd amphtude for the single soft-gluon emission (cf. (10)) 

r{rW') = /' -a,{xTWy{xTW' - Ql) , (37) 

71 T Jt X 



and A is the normahsation factor: 

r.l/3 



A(W-^g^) = exp 







dTV{TW 



2> 



(38) 



Qq is the lower bound of the transverse momentum of the emitted gluon, below which 
the soft gluon does not contribute to the particle proliferation. 

The multiplicity at given r is equal to the IPI cross section at given r integrated over 
the whole phase space and divided by dci/dr. The integrated IPl cross section at a given 
r (divided by a) is obtained by substituting the multiplicity Mg{xTW'^) for one of the Fs 
in each term of the series on the rhs of (36), thereby substituting l/(j — 1)! for the Bose 
factor in the j-th term: 

F(t) <n>r = — <n>r 
a dr 

= A{w\ Ql) [m'{tW^) + r ^M'(TiW^2)r ((1 - n)) w^) 
y Jo Ti 

+^ £ ^ dT2M'{nW')r{r2W')r ((i - n - r2)w') + • • j 
= 1^ ^2M'{nW^) {din - r) + r ((r - n)]^) 
+^iydT2r{T2W^)r{{i-n-T2)w^) + ---^ , (39) 



where 



M'{tW^) = =^ [ —as{xTW^)Mg{xTW^) . (40) 
Comparing the rhs of (39) with (36), we obtain 
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In (37) and (40), we included only the leading order term. The next-to- leading order 
corrections can be included by substituting 




xe{xTW^ -Ql 
^(2) ^C^(67_n^_5Nj_ 
2 V 6 2 9 



for (37), with two-loop a^, Aqcd = ^JTs'i 

M'{tW^) = ^ da; f- - 1 + ^) as{xTW^)Mg{xTW'') (43) 

for (40)0. 

The distribution -F(t) can be evaluated either directly (numerically) |]l4||jl5|, or using 
Laplace transformation [|l9l. It is an infrared-safe quantity in the following sense: Though 
each term of the series (36) depends on the infrared cutoff Qo, the dependence disappears 
in the sum at the high energy limit. In fact, the result somewhat depends on the value 
of Qq at a finite energy, which gives rise to the ambiguity of the order of Qi^jW . We 
do not know theoretically the exact value of Qo, because it is associated with the non- 
perturbative hadronisation process. 

In order to minimise the theoretical ambiguity, we use the experimentally observed 
thrust distribution for Fir) on the rhs of (41). Detailed data in the small r region was 
obtained by TOPAZ collaboration at TRISTAN at Vr=58GeV fg^. The result of the 



numerical evaluation of (41) is shown in Fig.3. Note that the experimental data of Fir) 
are given as the average values in finite bin sizes. Though the numerator of the rhs of 
(41) is continuous in r, its denominator is not. In order to reduce the error, we evaluate 
(41) only at the middle value of r in each bin, where the data of -F(t) is not far from its 
real value. In the integrand on the numerator. Fir — t\) is regarded as a step function 
(constant in each bin). In Fig.3, the result at each r is connected one another by straight 
lines. The curve thus obtained is still not smooth because of the inaccuracy included in 
this procedure. 

"^The term proportional to a^"^^ in the gluon emission amplitude off a quark (see (42)) gives the 
next-to-next correction in evaluating M'{tW'^). In fact, the leading-log approximation (36) is organised 
differently from MLLA, and we cannot tell a priori whether the next-to-next corrections in MLLA are 
in fact smaller than the next-to-leading order corrections in (36). (The former haven't been calculated 
yet.) 
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Let us now examine the validity of the approximation. First, we discuss the small r 
region. 

It has sometimes been argued that at the current energies the thrust distribution in the 
small T region (r < 0.07 or so) was determined mostly by the hadronisation process. The 
argument was partly based on the observation in the Monte Carlo simulation studies: The 
thrust distribution at the parton level (with partons produced above certain transverse 
momentum cutoff of the order of IGeV) is modified substantially in the small r region 
by the hadronisation (i.e. with the phenomenological hadronisation models to convert 
the particles produced by the QCD parton shower into hadrons). If so, use of (36) for 
F(r) (also (39) for (n)r) would not be justified in the small r region. 

A recent study , however, has shown that when we lower the cutoff in the partonic 
thrust distribution, the result reproduces the data quite well. It suggests either that the 
non-perturbative hadronisation process is almost irrelevant to the distribution, or that 
the non-perturbative hadropro duct ion can be simulated by the perturbative soft-gluon 
emission. Whichever the case, we may safely assume that the LPHD holds also in small 
r-processesQ. 

Another question to be addressed is on the recoil effect. In deriving the formulae, 
we have neglected the fact that an emitted gluon carries out a part of the energy of the 
parent quark, and accordingly the succeeding emission makes a reduced contribution to 
T (recoil effect) The recoil is large at the large r region. At the small r region, 

where soft-gluon emission dominates, the recoil is negligible. In fact, the correction at 
the large r region modifies the distribution at a small r through the change in its overall 
normalisation. The change in the normalisation of F{t), however, is canceled on the rhs 
of (41) between the numerator and the denominator. 

At large r, on the other hand, (41) overestimates the multiplicity. The thrust distri- 
bution in the large r-region is determined by the Feynmann diagrammes at the lowest 
orders in ag- Because the 0{as) matrix element vanishes for r > 1/3, the thrust dis- 
tribution decreases rapidly above r ^ 1/3. The numerator on rhs of (41), on the other 
hand, does not vanish at r > 1/3 even if F{ti) vanishes for ri > 1/3 in its integrand. 
The overestimation is particularly large when r is close to 1/3. 

Apart from this obvious overestimation near r = 1/3 (see Fig. 3 at r ^ 1/3), the 
recoil effect causes the reduction of the multiplicity in a wider r region. In addition to 

^ The assumption that the LPHD is correct with a small infrared cutoff is consistent with the expe- 
riences in the particle spectrum [Q. 
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the recoil effect on F{t), already discussed above, the recoil from the emission of a hard 
gluon reduces the multiplicity from the quark-antiquark pair, which partly compensates 
the multiplicity increase from the hard gluon (see the discussion given at the end of 
Set. 2). This effect is not included in (41). With the recoil effect, the multiplicity increase 
in r at large r would be less than the one estimated by (41). 

For the intermediate value of r (r ~0. 1-0.2), the overestimation reduces. On the 
other hand, we note that at very high energies the integrand in the numerator on the 
rhs of (41) has two sharp peaks, one at Ti (near the lower bound) and the other at 
Ti ~ T (near the upper bound). The contribution from the former peak corresponds to 
the multiplicity in the quark- (or antiquark-)jet, while the contribution from the latter 
corresponds to the multiplicity in the gluon-jet. We may therefore expect that (41) would 
give a prediction close to (20) in the intermediate r region. 

In Fig. 3, we also show {n)r evaluated by (20) in order to compare it with (41). We 
find that (41) gives smaller multiphcity than (20) (except near r = 1/3), despite the 
fact that the neglect of the recoil effect tends to give an overestimation in (41) in the 
large and intermediate r region. The reason for the difference is not difficult to find. 
At 58GeV, the thrust distribution is not very sharply peaked at very small r, so that 
in the integration in the numerator on the rhs of (41) the peak corresponding to the 
multiplicity in the gluon jet does not develop enough to give a contribution separate 
from the multiplicity from the qq-psdi. (The maximum of d(7/dln(l/T) is at r ~ 0.043.) 
If the jets are not collimated enough, it may be difficult to determine the multiplicity 
ratio from the multiplicity data with a thrust cut at this energy. 

Finally, let us make a comment on the determination of the normalisation of the 
multiplicity function Mg, which has been implicit so far. Because F{t) is a normalised 
distribution, we should have 



On the other hand, when we integrate F{T){n)r given by (41) over < r < 1, we obtain 




1/2 



(44) 



C drFir) 
Jo 



<n>T = 



1 



dr / 
Jo 



dri 



M'{tiW^)F{t - Ti) 



Jo Jo 



Ti 

/■I dri 



M'inW^) [ dT2F{T2)6{r - n - T2) 



Ti 



Jo 
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Jo Ti 
= < n >e+e- ■ 



(45) 



The final equality in (45) is correct in the accuracy of the next-to- leading order (see (34) in 
the last section). Indeed, the value of r cannot exceed 1/2, and therefore the condition 
(44) is not satisfied if the multiplicity function Mg is so normalised that the next-to- 
leading order formula for (n)e+e- (34) be satisfied. Alternatively, we may normalised 
the multiplicity function Mg so that the condition (44) be satisfied. We used the latter 
normalisation when we obtained the curve shown in Fig. 3. 

The condition (44), however, does not completely remove the error in the normalisa- 
tion. As is discussed above, the factorisation formula (41) gives overestimations of the 
multiplicity in the large r region. If the multiplicity is normalised by the use of (44), 
therefore, the multiplicity is underestimated in the small r region. When we compare the 
prediction with the experimental data in the small r region, the error in the normalisation 
(of the order of ag) should be taken into account. 

We thank M.Yamauchi at KEK for useful discussions. 
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W(GcV) 


r=0.1 


r=0.2 


r=0.3 


91 


1.33 


1.53 


1.60 


58 


1.31 


1.51 


1.58 



Table 1 A comparison of the predictions of the ratio (n)^/ (n)e+e- at different energies. 



FIGURE CAPTIONS 

Fig. 1 A schematic representation of the one-particle-inclusive cross section at a given r. 
The dashed hne represents the registered particle. Successive soft-gluon emission takes 
place in the shaded blob. 

(a) The registered particle originates from the hard quark. 

(b) The registered particle originates from the hard gluon. 

Fig. 2 The prediction for {{n)g + {n)qq)/ {n)e+e- (solid curve), and the multiplicity ratio 
R(t) = {n)g/{n)g (dashed curve). W = 58GeV and Aqcd = 0.15GeV. 

Fig. 3 The predictions for (n)^/(n)e+e- {W = 58GeV): The approximation (41) (sohd 
curve) vs the approximation (20) (dashed curve). 
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